Shear band formation in granular media as a variational problem 
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Strain in sheared dense granular material is often localized in a narrow region called shear band. 
Recent experiments in a modified Couette cell Q,Q provided localized shear flow in the bulk away 
from the confining walls. The non-trivial shape of the shear band was measured as the function of 
the cell geometry. First we present a geometric argument for narrow shear bands which connects 
the function of their surface position with the shape in the bulk. Assuming a simple dissipation 
mechanism we show that the principle of minimum dissipation of energy provides a good description 
of the shape function. Furthermore, we discuss the possibility and behavior of shear bands which 
are detached from the free surface and are entirely covered in the bulk. 
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Granular media constitute an interesting field of re- 
search from the point of view of both basic science and 
application. The intrinsically nonlinear and dissipative 
nature of the interaction between the particles leads to a 
great deal of interesting phenomena including force net- 
works, different kinds of instabilities, clustering, complex 
flow properties 0, 0] . One of the most apparent insta- 
bilities occurring in granular media is the formation of 
shear bands: At slow shear rate strain is not distributed 
throughout the sample but appears in a localized fashion 
along a rather narrow interface between two essentially 
unstrained parts. Shear-banding was the subject of var- 
ious experimental and theoretical studies in the last few 
years 

B H IS II @H E3 and stin presents significant 
difficulties for theoretical descriptions. 

Recently universal geometrical properties of shear 
bands were discovered in a Couette geometry modified 
such that shear localization near boundaries was avoided 
0,0- The experimental setup was a cylindrical container 
filled with grains up to height H. The bottom was split 
into an outer ring rotating with the container wall, and 
a stationary disk of radius R s in the center (Fig. [TJ 14] . 
Thus the outer and the inner part of the material were ro- 
tated relative to each other which created a shear band 
with cylindrical symmetry: It started at the perimeter 
of the stationary bottom disk and extended through the 
bulk up to the free surface. On the surface the angular 
velocity of the granular material as the function of the 
radius was measured. It follows to high accuracy an er- 
ror function characterized by two parameters: the width 
W and the center position R c of the shear zone. The 
width grows with increasing height H while the radius 
R c gets smaller. Interestingly, the surface position R c 
proved to be very robust: It depends only on two length 
parameters H and R s , but not on the particle properties 
nor on the shear rate. This is contrary to the width of 
the band, which is affected by the size and shape of the 
grains but is insensitive to changes of the slip radius R s . 



FIG. 1: Setup geometry. The rotating (white) and the sta- 
tionary (gray) parts induce shear flow in the granular material 
held by the container. 



The fact that R c and W depend on different control pa- 
rameters suggests that the two quantities can be studied 
separately. 

In this Letter we address the problem of the position 
of the shear band. By suitable choice e.g. of the particle 
size the width can be made arbitrarily small compared to 
H and R s . This justifies to model the shear band as an 
infinitely thin layer that represents the boundary between 
two blocks of material within which no flow occurs. We 
derive the shape of the band from a variational principle. 
Note that optimization principle was applied already to 
shear band formation in a different context |8j . 

Surface - bulk relation. In the experiments the po- 
sition R C (R S , H) on the free surface was found to scale 
like 



R c (R g ,H)=R s (l-(H/R a 



(1) 



with a Ri 2.5 for the experimentally accessible values of 
H/R s . It is much more difficult to measure the position 
of the shear band in the bulk, r(h), for fixed R s and H . 
Nonetheless, the experimental data clearly show that the 
bulk profile follows another form than R C (R S , H), and 
the bulk radius at height h depends also on the filling 
height H (Fig.©. 

We show that the bulk profile is determined once the 
surface positions R c {R Sl H) are given. Let us take a sys- 
tem with total height H and find the position of the shear 
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FIG. 2: The symbols are experimental data showing the 
shear zone radius r measured in the bulk at height h, taken 
from 0- +, x and o correspond filling heights 25, 37 and 
49 mm respectively, R s is 95 mm. The rectangle in the mid- 
dle shows the estimated errors in both directions for all data 
flU (plotted only for one data point). The solid line is the 
experimentally found fit curve for the surface positions. The 
dashed lines are the calculated bulk-positions based on the 
solid curve. 



band r at height h < H . The subsystem above h can be 
regarded as a smaller system with height (H — h) and 
with slip radius r at the bottom. Pressure and boundary 
conditions are the same, and the difference in the width 
is neglected in the narrow band approximation. We con- 
clude that 



R c (R s ,H)=R c (r,H-h). 



(2) 



Knowing the function R C (R S , H) thus allows to calculate 
the shape r(h) of the shear band throughout the whole 
system. Note, if the simple size scaling holds for the 
surface radius that R c /R s depends only on H/R s (as it 
was found in |2|) then it follows immediately also for the 
bulk function due to Eq. J2J: 
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(3) 



Based on Eq. an explicit functional form of the bulk 
profile can be obtained using the experimental fit func- 
tion of R r : 
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(4) 



The resulting curves for some filling heights are plotted 
in Fig. |21 using a — 2.5 and the comparison with the 
experimental data shows very good agreement. 

Variational principle. In order to describe the form 
of the shear band our idea is to apply the principle of 
least dissipation |12| (which is a common treatment of 
time-independent irreversible phenomena (l3j ). There- 
fore we require a steady state flow that matches the outer 
constraints but provides the minimum rate of energy dis- 
sipation. 

Applying this to the cylindrical geometry within the 
narrow band approximation the question of the shape is 



traced back to a variational problem among the functions 
r(h) (where H and R s are kept fixed) with the condition 
r(0) = R s while the other boundary at H is free. The 
dissipation rate is given by the sliding velocity between 
the two sides r(h)u> times the shear stress am integrated 
over the whole shear band. Up to a constant factor, the 
expression to be minimized is 



H 



r y 1 + (dr/dh) a tn dh 



(5) 



This quantity represents not only the dissipation rate but 
also the mechanical torque that the rotating and station- 
ary part of the system exert on each other. Therefore the 
least dissipation for this specific geometry is equivalent 
to the minimal torque which gives further justification 
of this approach: it is plausible that the yielding sur- 
face is established where the resistance against the outer 
constraint is the smallest, i.e. where the material is the 
weakest 

8]- 

Sliding model. For the shear stress in Eq. [51 we will 
use a very simple sliding model. The shear stress in the 
yielding surface is taken similar to the Coulomb friction 
between two solid bodies: It acts against the sliding di- 
rection, its magnitude is proportional to the normal pres- 
sure pressing the two sides against each other, but it is 
independent of the sliding velocity. We assume hydro- 
static pressure i.e. proportional to the depth ^5|. Thus 
Janssen effect is neglected, which is naturally justified if 
H is smaller than the container width. In our dynamical 
situation, however, we expect that the applicability of 
the hydrostatic pressure can be extended even for larger 
filling heights: The shear band (due to many collisions 
and slip events) acts as a source of small vibrations in the 
whole system and can cause slight creep at the particle- 
wall contacts inhibiting the particles to keep their original 
anchoring position. Finally they transmit their load to 
the next particle below rather than to the side wall and 
therefore the whole weight will be carried by the bottom. 

This sliding model leads to the variational problem 



H 



r yl + (dr/dh) (H - h)dh 



(6) 



The solutions which minimize the integral have auto- 
matically the scaling property given by Eq. © . The rea- 
son of this data collapse is that taking a A times larger 
system (i.e. taking Ai? s , XH, \r(h/X) instead of i? s , iJ, 
r(h)) changes the value of the integral only by a con- 
stant factor (A 4 ) thus it represents the same variational 
problem. 

Numerical results. The function r(h) is discretized 
and the minimization is performed numerically based on 
genetic optimization: r(h) is varied randomly but only 
the changes that lower the value of the left hand side in 
Eq. © ,i.e., reduce the dissipation, are admitted. During 
the optimization the noise level is continuously decreased 
and the final state r(h) is regarded as one local mini- 
mum of the variational problem. The landscape where 
the minimum has to be found is simple (see later). 




FIG. 3: Results obtained from the variational principle. Sym- 
bols show bulk profiles, from top to bottom H/R s — 0.15, 
0.25, 0.35, 0.45, 0.55 respectively. The two lines denote the 
surface positions as the function of the total height. The 
dashed line comes from our model, the solid line is the exper- 
imental fit function. 



Results shown in Fig. © reproduce nicely the qualita- 
tive behavior found in the experiment: the concave shear 
bands appear in the bulk and build up a convex confin- 
ing shape of the surface positions as the filling height is 
varied. The shear radii at a fixed bulk height h and also 
at the top get smaller with increasing filling height. 

The realistic bulk profiles provided by the principle of 
minimum dissipation can be interpreted easily. The bene- 
fit gained by having such a curved cylindrical shape (slim- 
mer at the top) in place of a regular cylinder is twofold: 
First, the surface of the shear band can be reduced by 
letting the regular cylinder deform and by pulling the 
shape towards the center a bit at a fixed bottom radius, 
similarly to a soap membrane spanned between a ring 
and a plate where the plane of the ring is parallel to the 
plate. Second, smaller radius results in smaller sliding ve- 
locity (or thinking of minimizing the torque it represents 
smaller lever). Therefore one can roughly think of the 
bulk profiles shown in Fig. © as equilibrium situations 
where the reduction of the sliding velocity is counterbal- 
anced by the increase in the surface due to going beyond 
the minimum surface. 

The quantitative agreement with the experimental fit 
function (Eq. ^) is also surprisingly good given the crude 
assumptions we made and the fact that our model con- 
tains no free fit parameters. The difference of the the- 
oretical and the experimental values of R c is less than 
20 % of R s - R c . 

We can more easily analyze the limit H/R s — > than 
in the experiment, where -Revalues are limited by the 
container size and a power law fit to the data at low 
filling heights is not very precise. From our variational 
principle we find an exponent a = 2 for small H/R s . In 
the region 0.4 < H/R s < 0.7 there is no clean power 
law any more and the effective exponent increases. In 
our model this seems to be due to a phase transition at 
about H/R s » 0.7 to a different shape of the shear band, 
which we discuss next. 

What is the predicted behavior of the shear band for 
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FIG. 4: Open and closed shear bands with cylindrical sym- 
metry, (a) upper radius, (b) bulk profiles for several filling 
heights (0 < H/R s < 7), (c,d) height of the shear band ht op - 



filling heights larger then those reported in 0, ? For 
large enough H/ R s the class of the "open" solutions dis- 
cussed so far is replaced by a new type of solutions: The 
shear band, instead of running up to the free surface and 
having a circle on the top as its upper edge, closes form- 
ing a cupola- like shape (still with bottom radius R s ). In 
that case the material covered by the "closed" shear band 
is at rest while the material around and above (includ- 
ing the whole free surface) is rotating. Several "open" 
and "closed" profiles can be seen in Fig. (0Jd) obtained 
for various values of H . Fig. (0J,) shows the upper radius 
of the shear bands which is characteristic of the "open" 
solutions but becomes zero for the "closed" ones. For 
these cupola shapes |l6j | a more relevant parameter is 
their heights h top < H in the center, plotted in Fig. (j3J;). 
For "open" profiles h top equals simply the system height. 

Interesting is the behavior of the parameter h top . It is 
a monotonically decreasing function of the filling height 
once it is detached from H, i.e. large filling heights 
press the shear band to the bottom. Solving the Euler- 
Lagrange equation for the variational problem neglecting 
terms of higher than first order in dh/dr and taking the 
limit of H 3> h t0 p one obtains 



h(r) = /Hop - 



Hop 



1 R 2 



(7) 



it is 



Fig. J3Ji) shows the numerical solution of htop 
in excellent agreement with the approximate analytical 
solution. In the limit H/R s — > oo the material forms one 
solid block and the sliding occurs on the surface of the 
bottom disk. 

The H dependence of htop has the following reason: 
Larger h corresponds to lower pressure and thus to 
smaller shear-resistance; therefore it is worth to raise the 
shear band into a cupola-shape even if its surface becomes 
larger than the disk at the bottom. Note, however, that 
stronger gravity would not affect the shape: It gives a 



4 





(b) - 













0.65 0.7 0.75 H 0.65 0.7 0.75 H 

FIG. 5: Hysteresis emerging in the variational problem. Up- 
per radius (a) and height (b) of the shear band in units of R B - 
The dashed line denotes the transition of the global minimum. 

larger pressure gradient but this constant factor has no 
impact on the variational problem. By contrast, increas- 
ing H (or, alternatively, applying additional pressure at 
the surface) makes the relative pressure change smaller 
near the bottom which results in weaker uplifting "force" . 
This is why the cupola height h top approaches zero for 
large H (respectively large pressure). 

At intermediate H (around 0.7R S ) there exists a small 
region of filling heights where the variational problem 
(Eq. H3) exhibits two local minima. The behavior of the 
local minima corresponds to a first order phase transi- 
tion, where the two phases are the two types of shear 
band profile (FigEJ. Following the "closed" solution as 
H is decreased its height grows till the top of the cupola 
reaches the surface at system height H-y. At this point 
the "closed" solution becomes unstable, R c at the top 
runs suddenly to a non-zero value (Fig. |5p,). Regard- 
ing the other direction (increasing H) the upper edge of 
the "open" shape is pulled towards the center then at 
height H2 it shrinks to one point which is followed by a 
jump into a "closed" shape (discontinuity in /i top see in 



Fig. 0}) . The interval [H\ , H2] defines the region where 
two local minima exist, outside this region there is only 
one phase possible. Local minima are physically mean- 
ingful in the presence of a kinetic barrier. The latter 
can be provided by the difference between the static and 
dynamic friction coefficients, usually present in granular 
systems. 

In this Letter we have presented a theoretical analy- 
sis of recent experiments [1,2] on shear band formation 
in granular media. We used the approximation of nar- 
row shear bands. First we showed a geometric argument 
which related the position of the surface endpoints of 
the band to the bulk shape. We calculated the shape of 
shear bands from a variational principle and the results 
are in good agreement with the experiments. The theory 
provides with a number of predictions, i) There is a tran- 
sition in the shape of the shear band as a function of the 
filling height H : For low values of H the shear band is an 
open curved cylinder which ends on the surface while for 
large H a cupola is formed, ii) This transition is of the 
first order, accompanied by a hysteresis, iii) The height 
of the cupola is proportional to R s /H . These predictions 
should be experimentally accessible. 

An interesting question we did not address in this Let- 
ter is the nuclcation kinetics at the first order phase tran- 
sition from open to closed shapes. It is conceivable that 
the finite width of the shear bands plays a role here. 
Further progress in this direction requires a continuum 
theory going beyond the narrow band approximation we 
used here. 
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